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, Abstract. We use Fabry-Perot Ha spectroscopy, complemented with 

published H I radio synthesis observations to derive high resolution rota- 
, tion curves of spiral galaxies. We investigate precisely their inner mass 

' distribution and compare it to CDM simulations predictions. Having 

. verified the existence of the so-called core-cusp problem, we find that 

the dark halo density inner slope is related to the galaxy masses. Dwarf 
'q_|. galaxies with Vmax < 100 km/s have halo density inner slope < 7 < 0.7 

^ [ while galaxies with Vmax > 100 km/s are best fitted by 7 > 1. 

u 
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1. Introduction 



The dynamical masses of spiral galaxies are known to differ significantly from 
\ their visible masses. The commonly accepted cause is the existence of a ellip- 

soidal halo of unseen matter in addition to the stars and gas. The exact density 
distribution of these halos have become an increasingly important issue. On one 
side, N-body simulations of the cosmological evolution of the Cold Dark Matter 
(CDM) have now reached a sufficient resolution to predict the dark halo density 
profiles down to a scale corresponding to the inner parts of the spiral galaxies 
(Fukushige et al. 1997; Navarro et al. 1997; Moore et al. 1998; Ghigna et al. 
2000). They almost always predict dense cuspy halos. On the other side, obser- 
vations of dwarf spiral galaxies show shallow inner rotation curves, compatible 
with a fiat density core (Blais-Ouellette et al. 2000, hereafter paper H). 

The first step in showing the reality of this discrepancy is to eliminate the 
known possibilities of systematic observational biases. Two classes of errors 
could contribute to underestimate the velocities, hence the computed density, 
in the inner parts of spiral galaxies. The prime culprit in radio observations 
is the "beam smearing" effect due to the relatively low angular resolution of 
21 cm data with sufficient sensitivity to detect H I in the outer part of spiral 
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galaxies. Combining the H I density gradient with the velocity gradient will lead 
to underestimate the velocity at a given radius. 

Ha observations always easily reach an angular resolution where any beam 
smearing effect can be neglected. A less often commented source of uncertainties 
though is found in long slit observations, where most of Ha data come from. The 
lack of 2D coverage makes the alignment of the slit crucial to retrieve the real 
kinematics of a galaxy. Missing the kincmatical center, which is not always 
the photometric center (paper II), or just a few degrees between the slit and 
the galaxy position angle will also lead to an underestimation of the velocities. 
Inclination estimation, which has to be photometrically determined, is another 
major source of uncertainties. In addition, the presence of a bar would hardly 
be noticed and its effect would most probably be confounded with the rotational 
kinematics. 

In order to hedge ourselves against these biases, wc use Fabry-Pcrot high 
resolution Ha spectroscopy combined with published radio synthesis observa- 
tions to study the mass distribution of 8 dwarf and spiral galaxies. In paper 
II, with a smaller sample, wc focused on modeling the mass distribution using 
different shape of dark matter halos. Here, in addition, we address the more 
precise question of what inner density slope (7) halo can have for a given galaxy 
type or mass. 

In section 2, we first briefly review the mass modeling used in the study. 
Then, in section 3, we look at detailed mass models of a few galaxies. In section 
4, the relation between 7 and galaxy mass is discussed followed by concluding 
remarks. 

2. Modeling the Mass Distribution 

To investigate in details the mass distribution of dark matter halos without a few 
assumptions on the matter content of spiral galaxies, one would have to adjust 
a good dozen of parameters. First, the luminous matter distribution depends 
on the disk and bulge mass-to-light ratios (and their radial gradient), and on 
the bulge-to-disk ratio. The H I contribution have to be corrected for helium 
fraction. The dark halos can be non-spherical, in addition to the five parameters 
usually used to describe the radial density distribution. This general distribution 
function can be expressed as (see Paper II for details): 



where po ^'^'^ '"o ^"^^ characteristic density and radius, where c, included for ease 
of comparison with other works, can force the presence of a flat density core, 
where a and 7 are respectively the inverse outer and inner logarithmic slopes, 
and /3 the transition parameter. 

One could add the distance that is used to calculate the light distribution, 
and a central mass which is suspected to exist in most spirals. 

For some of these constraints, a fixed value is well accepted. The mild or 
absent color gradient in spiral galaxies lead to a radially constant mass-to-light 
ratio. The helium fraction can be approximate to its primordial abundance. 
Distances are hopefully well constrained by independent means. 
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Otherwise, essentially three data sets can be used to constraint these param- 
eters: luminosity distributions (visible an H I ) and velocity field or, assuming 
axisymmetry, light profiles and rotation curve. The visible light profile, is used 
to determine the bulge to disk ratio. The other parameters are all to the charge 
of the rotation curve. Most of the luminous contributions are heavily constrained 
by the most inner parts of the curve leading to a possible degeneracy between 
central mass, bulge, inner disk and halo contributions. That is why, many stud- 
ies including this one, tend to focus on dwarf galaxies where bulges are negligible 
and large central masses excluded by the rotation curve. Only there, can one 
put strong limits on the inner slope of dark halo density distribution. For these 
galaxies, luminous matter is dependent only on the disk mass-to-light ratio, and 
dark matter from the central density, the core radius, and the three shape pa- 
rameters. From the latter, only 7, the inner slope, has a significant impact on 
halo shape at the scale of dwarf galaxies, a and /3 being very poorly constrained. 

In our procedure, the visible light profile (eventually decomposed in its bulge 
and disk components) is inverted in a mass profile scaled by the mass-to-light 
ratios. The H I contribution is estimated from its distribution and multiply by 
1.33 to account for helium. The halo shape is, in general, fixed to a known 
profile. We choose four profiles, widely used in the literature, that span, with 
some redundancy, through most possibilities. The cuspy profile from Navarro 
et al. (1997) is a canonical example of CDM simulation prediction, the pseudo- 
isothermal sphere, the profile from Kravtsov et al. (1998), and Burkert's profile 
(Burkert et al. 1995) are all cusp-less profile with different shape parameters. 
Rotational velocities from all the contributions are then added in quadrature 
and compare to the rotation curve using a standard minimization. 

3. Mass Models 

3.1. The Sample 

Beside the dwarf galaxies NGC 3109 and IC 2574 from paper II, the whole 
sample includes 3 fairly bulge-less spirals (UGC 2259, NGC 2403, and NGC 
6946), and 3 earlier type spirals (NGC 5055, NGC 2841, and NGC 5985) from 
paper III. Details are given in Table |l|. 

3.2. The Results 

Complete modeling, using the four models, of the most relevant cases are pre- 
sented here. NGC 3109 and IC 2574 (from paper II) are re-analyze to correct 
for a numerical problem which underestimated the central density. They are 
good examples of CDM incompatible dwarfs. Are also modeled, NGC 2403, 
as an example of bulge-less small spiral compatible with any model, and NGC 
5055, as an earlier type spiral with a bulge. Analysis of the whole sample will 
be presented in a future paper (paper IV). Error bars are the quadratic sum 
of: half the velocity difference between receding and approaching sides, half the 
correction for asymmetric drift, error cause by uncertainty on inclination, and 
statistical error (cr/\/iV). 

It is rather clear that NFW profile is incompatible with the two dwarf 
galaxies. These results are in line with most similar studies using H I (e.g. Cote 
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Figure 1. Best fit mass model for NGC 3109 (using the Ha up to 
410"completed by the H I ), and for IC 2574 (H I only). 
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Figure 2. Same as Fig. | for NGC 2403 (using the Ha up to 345"), 
and for NGC 5055 (Ha is used up to 52 "). 
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Table 1. Parameters of the sample. From the main references unless 
otherwise specified. 



Name 


Type 


Distance 


D25 




a-' 


Ms 


Li3 




References 






Mpc 






11 




IO8L0 


H I 


Ha 


photom. 


IC 2574 


SABm 


3.0 


9.76 


8.63 


165 


-16.87 




m 


cc 


m 


NGC 3109 


SBm 


1.36" 


14.4 


13.3 


234.9 









cc 





NGC 5585 


SABd 


6.2 


5.27 


3.62 


46.7 






P 


gg 


P 


UGC 2259 


SBcd 


9.6'= 


2.6 


1.9 


28.2 


-16.33 


5.06 


b 


ii 


b 


NGC 2403 


SABcd 


3.2'' 


21.88 


13.0 


134 


-19.50 


7.9 


c 


ii 


1 


NGC 6946 


SABcd 


7 


11.5 


7.8 


1.92 


-21.38 


530 


d 


ii 


d 


NGC 3198 


SBc 


9.15 












c 


nn 


1 


NGC 5055 


SAbc 


9.2'= 


12.56 


9.8 


108.9 


-21.15 




e 


ii 


1 



"This study 
"'Martimbeau et al. 
"Musella et al.97 
'^'^paperll 
°Jobin et al. 1990 
PCote et al. 1991 



1994 



1988 



^^paperl 
''Carignan et al. 
"paperlll 

''Freedman et al. 1990 
^Begeman et al.l987 
'Kent et al. 1987 



"Carignan et al. 1990 
""Corradi et al. 1991 
'^Thornley et al. 1997 
*Aaronson et al. 1983 
'=Ho = 75kms-VMpc 



et al. 2000) or long slit observations (e.g. Swaters et al. 2000, de Blok et al. 
2001) but without the related uncertainties (van den Bosch et al. 2000). It 
has to be noticed that higher resolution N-body simulations tend to give even 
steeper inner density slopes. The present discrepancy is therefore genuine and 
probably involve subtle phenomena or new physics. 

Mass-to-light ratios tend to be unrealistically low for the smallest galax- 
ies. Halo contributions are therefore an upper limit and are in fact probably 
shallower. 



4. Halo Density Gradient 

The range of possible observational biases that could explain the core-cusp prob- 
lem is now significantly reduced, rather close to nil. The question of the physical 
cause of this shallow inner density distribution have been addressed many times 
leading to proposals ranging from multi-component dark matter (Burkert & Silk 
1997) to self-interacting dark matter (Spergel & Steinhardt) and baryonic feed- 
back processes. These scenarios are though to behave differently under different 
gravitational potentials. For example, feedback processes can hardly have an 
impact in galaxies with very deep potential. Plotting 7, the density profile inner 
slope, against the asymptotic rotational velocity would precisely show the evo- 
lution of this behavior. Figure ^ plots 7 against the observe maximum velocity 
which, for NGC 3109 and IC 2574, is lower than the true asymptotic velocity. 

The steepness of NGC 5055 rotation curve leads to an important degeneracy 
between 7, ro and rho^. It is therefore not clear from the plot if there is a canonic 
value of 7 around 1.2 that breaks down below 100 km/s or if 7 is intrinsically 
increasing in massive galaxies. 
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Table 2. Parameters of the mass models 



Model 


Galaxy 


Type 


(7W/Lb)* bulge 


(7W/Lb)* disk 


ro 
kpc 


Po 

10 M.Q /pc 


'J 

x" 


ISO 


NGC 3109 


SBm 


no bulge 


0.1 


2.4 


2.4 


0.44 




IC 2574 


SABm 


no bulge 


0.3 


5.0 


0.75 


2.7 




NGC 2403 


SABcd 


no bulge 


2.5 


4.8 


1.7 


1.9 




NGC 5055 


SABbc 


2.3 


0.8 


8.2 


3.4 


0.55 


Burkert 


NGC 3109 




no bulge 


0.4 


3.0 


1.8 


0.22 




IC 2574 




no bulge 


0.3 


8.9 


0.8 


2.5 




NGC 2403 




no bulge 


2.5 


8.0 


1.9 


1.8 




NGC 5055 




3.2 


1.6 


8.2 


2.8 


.56 


KKBP 


NGC 3109 




no bulge 


0.0 


4.0 


1.7 


0.5 




IC 2574 




no bulge 


0.3 


10.5 


0.46 


2.5 




NGC 2403 




no bulge 


2.5 


8.5 


1.2 


1.8 




NGC 5055 




3.1 


2.0 


11.2 


1.0 


.56 


NFW 


NGC 3109 




no bulge 


0.0 


605 


0.0035 


9.0 




IC 2574 




no bulge 


0.0 


295 


0.005 


32 




NGC 2403 




no bulge 


1.4 


11.8 


0.9 


1.4 




NGC 5055 




2.6 


1.1 


10.8 


1.8 


.52 



ro : core radius of the dark halo 
Pq: central density of the dark halo 



5. Conclusion 

We have modeled 8 dwarf and spiral galaxies using cusp-less and cuspy halos. 
The latter, predicted by CDM N-body simulations are clearly incompatible with 
dwarf galaxy kinematics. More precisely, galaxies with less than 100 km/s of 
maximum rotation velocity have a inner density logarithmic slope (7) of less 
than .7 as opposed to 7 > 1 predicted by the N-body simulations. 

Presently, a, the outer density slope is poorly constrained by rotation 
curves. One need to reach a radius where the luminous disk contribution is 
negligible while the rotation curve is fiat or decreasing (Carignan h Purton 
1998). In our sample, UGC 2259 is a good candidate. 

Future studies should in part extend to massive bulge-less galaxies to see if 
7 is intrinsically rising with galaxy mass or if it is constant, with a break down 
at low mass. The GHASP survey (Garrido &; Amram 2002) should allow this 
kind of study. 
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Figure 3. Inner logarithmic slope (7) of the CDM halos. Note that 
N-body simulations give 7 > 1. 



Carignan, C. & Freeman, K.C. 1988, ApJ, 332, L33 

Carignan, C, Charbonncau, P., Boulangcr, F. & Viallefond, F. 1990, A&A, 234, 43 

Carignan, C. & Purton, C. 1998, ApJ, 506, 125 

Cote, S., Carignan, C, h Freeman, K. C. 2000, AJ, 120, 3027 

Cote, S., Carignan, C. & Sancisi, R. 1991, AJ, 102, 904 

Corradi, R.L.M., Boulesteix, J., Bosma, A., Capaccioli, M., Amram, P. & Marcelin, M. 
1991, A&A, 244, 27 

de Blok, W. J. G., McGaugh, S. S., Bosma, A. & Rubin, V. C. 2001, ApJ, 552, L23 
Freedman, W. L. 1990, ApJ, 335, L35 
Fukushige & T., Makino, J. 1997, ApJ, 477, L9 
Garrido, O., Amram, P., Tiiis conference 

Ghigna, S., Moore, B., Governato, F. ,Lake, G. ,Quinn, T. & Stadel, J. 2000, 544, 616 
Jobin, M. & Carignan, C. 1990, AJ, 100, 648 
Kent, S. M. 1984, ApJS, 56, 105 

Kravtsov, A. V., Klypin, A. A., Bullock, J. S. & Primack, J. R. 2000,ApJ, 502, 48 

Martimbeau, N., Carignan, C. & Roy J.-R. 1994, AJ, 107, 543 

Moore, B., Governato, F., Quinn, T., Stadcl, J. & Lake, G. 1998, ApJ, 498, L5 

Musella, I., Piotto, G. & Capaccioli, M. 1997, AJ, 114, 976 

Navarro, J.F., Frenk C.S. & White, S.D.M. 1997, ApJ, 490, 493 

Swatcrs, R. A., Madore, B. F. & Trewhella, M. 2000, ApJ, 531, L107 

Thornley M. D. & Mundy L. G. 1997, ApJ, 484, 202 

van den Bosch, F. C, Robertson, B. E., Dalcanton, J. J. & de Blok, W. J. G. 2000, AJ, 
119, 1579 



